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On what scale should inflationary observables be constrained? 
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We examine tlie clioice of scale at wliicli constraints on inflationary observables are presented. 
We describe an implementation of the hierarchy of inflationary consistency equations which ensures 
that they remain enforced on different scales, and then seek to optimize the scale for presentation of 
constraints on marginalized inflationary parameters from WMAP3 data. For models with spectral 
index running, we find a strong variation of the constraints through the range of observational scales 
available, and optimize by finding the scale which decorrelates constraints on the spectral index ns 
and the running. This scale is fe = 0.017 Mpc~^, and gives a reduction by a factor of more than 
four in the allowed parameter area in the ns-r plane (r being the tensor-to-scalar ratio) relative to 
k = 0.002 Mpc~^. These optimized constraints are similar to those obtained in the no-running case. 
We also extend the analysis to a larger compilation of data, finding essentially the same conclusions. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

In presenting observational constraints on the primor- 
dial power spectra, such as those that may have been 
generated by inflation, it is necessary to specify the scale 
at which the observables are being determined. Provided 
the full posterior distribution over all parameters is given, 
this choice is an arbitrary one. However, if the informa- 
tion is to be compressed via marginalization, the choice 
of this scale matters, and should be chosen in order to 
optimize the presentation of constraints. 

In the Wilkinson Microwave Anisotropy Probe 
(WMAP) three-year cosmological parameters paper [l| 
the scale 0.002 Mpc~^ is used, which is close to the ob- 
servable horizon, while Kurki-Suonio et al. f2j and Finelli 
et al. 3 found that the choice of 0.01 Mpc~^ worked bet- 
ter in constraining inflationary observables as it is closer 
to the statistical center of the data. The scale 0.05 Mpc^^ 
is also commonly used, being the default scale of the 
CosmoMC package [l]- The pivot scale was also dis- 
cussed in Ref. @, who sought the scale where the per- 
turbation amplitude was best determined (decorrelated 
with other power spectrum parameters), and in Ref. 
who sought the scale at which the perturbation spec- 
trum reconstructed using the flow formalism was best 
constrained. 

In this paper we make a systematic exploration of the 
choice of scale in the context of inflation models. This 
choice is particularly important in cases where models 
with significant spectral index running are allowed. Such 
models have received quite a bit of attention since the 
WMAP results emerged (see e.g. Refs. 0,[l[i|)- 



II. METHODOLOGY 

For definiteness we concentrate on single-field infla- 
tionary models, though many of the issues we discuss 



are more general. These models predict spectra of scalar 
and tensor perturbations which are related by a hierar- 
chy of consistency equations (T^jIHIj the first of which is, 
at lowest-order, the well-known relation r = — Snx where 
r is the tensor-to-scalar ratio and tit the tensor spectral 
index. These parameters can in turn be related to the 
inflationary slow-roll parameters describing the shape of 
the potential. 

Our main aim in this paper is to examine the optimal 
choice of scale at which to present observational con- 
straints on inflation. In order to fit the spectra from 
data, they must first be parametrized, which is usually 
done by specifying their amplitude and some number of 
derivatives (i.e. the spectral index, running, etc) at a par- 
ticular scale. So far, this scale has been chosen by hand. 

The choice of scale, being arbitrary, ought not to affect 
the conclusions one draws. However this is only the case 
if one specifies the full multi-dimensional posterior pa- 
rameter distributions, and provided the model definition 
is internally self-consistent. The first of these is often 
not the case, as one commonly wishes to condense infor- 
mation into one- or two-dimensional marginalized pos- 
terior distributions, which throws away the information 
on parameter correlations necessary to translate between 
scales. The second condition of model self-consistency 
holds in most circumstances, but often not in the way 
infiationary spectra are implemented. 

The problem of model definition in inflationary models 
is the enforcement of the consistency equations between 
scalars and tensors. Typically both spectra are allowed 
to be power-laws but with different indices; if the usual 
consistency equation is enforced at one scale, it will then 
no longer hold at any other. Put another way, if the 
scalars are a perfect power-law, then the tensor spectrum 
implied by the consistency relations is not (unless the 
spectral indices are the same) . Yet another way, the set of 
models generated by imposing the consistency equation 
at one scale is a different set of models from that obtained 
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using another scale. This problem is further exacerbated 
if authors go on to include scalar spectral index running, 
while perhaps still leaving the tensors as a power-law. 
Before discussing the choice of scale, we should therefore 
first fix this problem (while admitting that the difference 
may be too small to be very important). 

This is achieved by implementing the full inflation- 
ary consistency equation hierarchy, as given explicitly in 
Ref. [ll|. As well as the first consistency equation, this 
enforces that each derivative of the consistency equation 
also holds at a given scale. When using a Taylor ex- 
pansion to shift from one scale to another, this hierarchy 
then ensures that the consistency equations will still hold 
at the new scale (up to some level set by the truncation 
of the hierarchy). 

We note that these complications are needed only if one 
fits the phenomenological parameters (amplitude, spec- 
tral index, running, etc) from the data and then trans- 
lates to infiationary observables. If instead one fits the 
slow-roll parameters directly (e.g. Ref. [I,[l3|) or via flow 
equations d, [l^ then the consistency equation hierar- 
chy is automatically enforced. 

We consider a parametrization of the scalar and tensor 
perturbations as follows 

Al{k) oc (fc/A:^)("s-i)+(<i«s/di„fe) i„fc/fc. 

A^{k) cx (fc//c^)"T+(<i«T/di„fe) i„fc/fc. ^ (2) 

the constants of proportionality being the amplitude of 
the perturbations at scale fc*. The tensor-to-scalar ra- 
tio is deflned by r{k) = 16A^(fc)/Ag(fc), and the ten- 
sor spectral index is determined via the flrst consistency 
equation. 

In order that the first consistency equation is enforced 
at all scales (to linear order in A In A;), we need to imple- 
ment the second consistency equation to fix the tensor 
running, which is not a genuine new degree of freedom. 
This second equation is given by ITll | 



driT 
dink 



= TiT [riT - {ng - 1)] 



(3) 



We enforce this when carrying out our data-fitting. 

One could further enforce higher consistency equa- 
tions, so that for instance the second consistency equa- 
tion also is preserved under change of scales. However 
current data quality is a long way from the point where 
doing so would make any practical difference, since the 
tensors are potentially observable only over a limited 
range of scales. 

We use the Monte Carlo Markov Chain (MCMC) tech- 
nique to explore the parameter space, using the Cos- 
moMC package W\ . We consider a ACDM model in a flat 
universe and take fc* = 0.05 Mpc~^ as the scale where all 
power spectrum parameters are defined when fitting to 
data. We vary up to eight parameters 



dns 



where fib^^ and fl^^h'^ are the physical baryon and dark 
matter densities, 9 is the ratio of the sound horizon to 
the angular-diameter distance, t is the optical depth, and 
the remaining parameters specify the power spectra. We 
apply a set of uniform priors: 



0.005 < n^h^ < 0.1 
0.5 < 61 < 10 
0.5 < ns < 1.5 

< r < 2 



0.01 < ildnih^ < 0.99 

0.01 < r < 0.8 
2.7 < log(10i''v4s) < 4 
-0.2 < dns/d\nk < 0.2 



Until Section |Vl our constraints are from WMAP3 data 
alone. 



III. CHOICE OF SCALE: MODELS WITH 
SCALAR RUNNING 

We first consider models which allow running of the 
scalar spectral index, which we will see is the case where 
the choice of scale is most important. For comparison, 
models without running are studied in the next section. 



A. Tilt and running 

The simplest combination of observables to consider is 
the tilt and running of the scalars. Observational im- 
plications of this were first discussed in Ref. [l3|, which 
forecasted CMB constraints from the Planck satellite on 
running spectral index models. The paper pointed out 
that there would be a scale at which the uncertainties on 
tilt and running would become uncorrelated, and that 
(at least in a gaussian approximation) on that scale the 
uncertainty in n would recover its value for the case of 
no running.^ This could be spoiled by degeneracies with 
oth er p arameters, but at Planck accuracy appears not to 
be 111. 

Anyway, we wish to find the scale at which the tilt 
and running decorrelate for actual current data. To do 
this we take the distribution of these two variables as 
given by the MCMC analysis, which specifies quantities 
at k^ = 0.05 Mpc~^. We then fit the chain elements with 
a linear relation, ng = A + B dns/dlnk, and by inserting 
into the expression 



ns{k) = ns{k^) 



dns J k 
dink k^ 



(4) 



we arrive at a condition for the difference in scale which 
decorrelates ng and dng/dlnk: B — — Infc/fc*. This 
scale turns out io he k — 0.017 Mpc^^. Then we use 
Eq. (HI) to convert the distribution at scale fc, to the one 
at scale k — 0.017 Mpc~^ to obtain the decorrelated ns 



dh\k 



^ This observation was actually credited to Daniel Eisenstein, who 
was not an author of that paper. 
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FIG. 1: Constraints in the ns-a plane (where a = dns/d\nk) 
at several scales, k — 0.017 Mpc~^ is the decorrelation scale 
for these parameters. 



and dns/dlnk. More generally, we can explore the con- 
straints at other scales via the same formalism. The con- 
straints at a set of different scales, including the WMAP 
scale and the decorrelation scale, are shown in Fig. [TJ 

In this particular case (but not those that follow), the 
transformation between parameters induced by the scale 
change has unit Jacobian, which means that the 2D con- 
tour areas are preserved. However it is clear from Fig. [1] 
that the marginalized uncertainty on ng at the decorre- 
lation scale is significantly smaller. The WMAP choice, 
0.002 Mpc~^, gives a significantly-angled constraint area 
and is clearly to be avoided if useful marginalized con- 
straints on ns are to be quoted. Unfortunately, the main 
WMAP3 results for models with running are presented 
at this scale. 

For our choice of parameters and dataset (WMAP3 
alone), from separate fits where running is not included 
we find the marginalized constraint on ng is ng — 



0.993 



+0.029,-1-0.067 
-0.030,-0.053 



(at 68% and 95% confidence). With 
running, the marginalized constraint at the decorrelation 
scale is ng = 0-981!^o o34'lo'o53. As anticipated, there- 
fore, when including running the shift in the best-fit ng at 
the optimized scale is negligible within the uncertainty. 
This is somewhat trivial as it could have been chosen 
to match exactly by specific choice of scale — choos- 
ing k = 0.015 Mpc~^ achieves this. Much more impor- 
tantly, we see that the uncertainty on ng at the decorre- 
lation scale is hardly increased when running is included, 
whereas it is greatly increased at e.g. 0.002 Mpc~^. The 
ID marginalized constraints on all parameters have min- 
imum uncertainty at the decorrelation scale. 

Incidentally, for the scalar running the marginal- 



FIG. 2: Constraints on ns versus r at several scales. 



ized constraints we obtained are dng/dlnk 
-0.075 



-1-0.041,-1-0,082 
-0,043,-0,093' 



very similar to those quoted by 



WMAP3 for models with running and tensors 



B. Tilt and the tensor— scalar ratio 

We now turn to other combinations of observables, rel- 
evant to constraining inflation. 

To obtain r at other scales we perform an expansion, to 
the order considered, of the scalar and tensor amplitudes. 
The relation is 

r(k) l + »TlnA + l[4 + -^]ln^£ 

r{h) 1 + (ng - 1) In A + 1 [(^^ i)2 + ^] k_ ■ 

(5) 

where all observables without an argument '(k)' are eval- 
uated at = 0.05 Mpc~^, and where dnx/rflnfe was set 
according to the lowest-order version of the second con- 
sistency equation, Eq. ([3]). Having expressions for ng and 
r at different scales, we can now choose several scales and 
get the distribution of the two variables at each, shown 
in Fig. H 

In this case the transformation alters the contour areas 
as well as distorting them. The middle panel of Fig. [3] 
shows the areas enclosed by the 95% confidence contour 
in the ng-r plane at different scales.^ The top panel 



^ These values were obtained by taking the number of points in a 
50 X 50 grid that lie within that contour. The number of grid 
points across each axis corresponds also to the number of bins 
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FIG. 3: Variation of parameter plane area with scale. For 
the e-rj case both lowest (black, full) and next order (red, 
dashed) are shown. For ns versus running the area should 
be independent of scale, and the variations indicate the noise 
level in the area estimation. 



shows the same for ns and running discussed in the pre- 
vious subsection. In the ns-r plane the minimum area 
was near k = 0.017 Mpc ^ as expected (the precise value 
found was slightly smaller). As inflation model builders 
typically just look at these marginalized plots to decide if 
their model is viable, it is clearly important to present the 
constraints at a good scale. 0.002 Mpc~^ is not a good 
scale for this purpose, as has previously been stressed 
also in Ref. (al. 



C. Inflationary slow-roll parameters: lowest order 
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FIG. 4: Constraints on e versus rj, at lowest order, evaluated 
at several scales. 



The pivot scale fc* corresponds to some particular scalar 
field value (defined as the field value when fc* = aH 
during inflation) , in the vicinity of which the scalar field 
potential is being reconstructed. Shifting the pivot scale 
means expanding about a different point on the potential. 

We first concentrate on the constraints given at lowest 
order, taking the expressions for the potential at this 
order by Lidsey et al. [lol |: 
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(7) 



4(fc) 



3r 

2^' 



nsik)] 



(where without loss of generality we take (j) to increase in 
time). From these the first two slow- roll parameters are 
expressed in terms of the observables, to lowest order, by 



We now examine how the constraints on the first two 
slow-roll parameters e and 77 are affected by scale change. 
We take the usual definitions in terms of the potential |15| 



'"PI 
167r 



used to sample the distribution. We found that accurate area 
estimation needed at least 50 bins, though such an aggressive 
binning level leads to less smooth contours than are usually seen. 



16 ' 



3 1, 

— r (1 — riR) 

16 2^ ^' 



(8) 



Shifting the scale of the observables shifts the location 
on the potential, and at lowest-order the constraints on 
e and 77 then become independent of the running at that 
scale (which could be used to determine a third slow-roll 
parameter ^ = ml^/8TTy^V'V"' /V^). 

The results are presented in Fig HI and again show 
strong variation of the allowed parameter area with 
choice of scale, as indicated in Fig. [31 



5 




k=0.002 



0.02 0.04 



k=0.005 



0.02 0.04 




k=0.009 



0.02 0.04 



1 

0.8 
•- 0.6 

0.4 

0.2 



k=0.002 



0.8 1 1.2 1.4 



k=0.005 



0.8 1 1.2 1.4 



k=0.009 



0.8 1 1.2 1.4 



k=0.013 



k=0.017 





k=0.019 




0.02 0.04 



O.02 0.04 



0.02 0.04 



k=0.013 



0.8 1 1.2 1.4 



1 k=0.017 



0.8 1 1.2 1.4 



k=0.019 



0.8 1 1.2 1.4 



k=0.03 




k=0.05 



k=0.07 




0.02 0.04 



0.02 0.04 



0.02 0.04 



0.8 
•- 0.6 

0.4 

0.2 



k=0.03 



0.8 1 1.2 1.4 



k=0.05 



0.4 
0.2 



0.8 1 1.2 1.4 



k=0.07 



0.8 1 1.2 1.4 



FIG. 5: Constraints on e versus t], to next order, at several 
scales. 



FIG. 6: Constraints on ns and r when no scalar running is 
present. 



D. Inflationary slow-roll parameters: next order 

Now we can take the expressions for the potential to 
next order also given by Lidsey et al. [lO |: 



V'{cp) 



TSVtt 3 A\ 



'PI 



257r 



As 
' Al 



(36C + 2)^-i(l-nsf 



1 - 0.85-^3: + 0.53(1 - ?is) 
3 
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(9) 



(12C- 6)^(1 -ns)--(3C-l) 



din A: 



where C = -2 + hi2 + 7 ~ -0.73, 7 is the Euler- 
Mascheroni constant, and again the </) value corresponds 
to horizon crossing of the scale at which the constraints 
are being imposed. 

With these next-order expressions for the potential, e 
and r\ are 



r l-0.85r/16 + 0.53(1 -ns) 
16 1 + 0.21r/16 

1 1 r 9 3 

■ — { — r (1 — nc) 

3 l + 0.21r/16ll6 2^ ^' 

+ (36C + 2) (^)'-i(l-nsf 
-(12C -6)^(1 -ns)-^(3C - 



(10) 
(11) 



dln/c. 



The second parameter now depends on the running. The 
running term has a coefficient of about one half, and 
given how weakly running is constrained this term has a 
significant impact on the constraints. 

The constraints at each scale are presented in Fig [S] 
The picture here is rather different, with the area chang- 
ing much more slowly as h is decreased, and the minimum 
area being at a much smaller h. This is because for typ- 
ical models the next-order correction from the running 
happens to be comparable to the change in the lowest- 
order expression for 77 coming from the changing ns , also 
induced by the running, as the scale changes. These 
terms approximately cancel going to smaller fc, i.e. the 
constraints change less when simultaneously reducing fc 
and introducing next-order corrections than they would 
if only one of these were done. This is just a coincidence 
(and not much of a coincidence at that, since partial can- 
cellation would have to happen as k was changed in one 
or other direction) of no great significance, and will go 
away when in future running is better constrained. 



IV. CHOICE OF SCALE: MODELS WITH NO 
SCALAR RUNNING 

For comparison, we now take a look at models where 
no running of the scalar index is allowed.^ In this case 



^ We still keep the tensor running in the analysis, however. It is 
not an additional degree of freedom, its inclusion ensuring the 
validity of the first consistency equation at all scales. 
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FIG. 7: Comparison of constraints in the n^-r plane at the 
optimal scale with no running (dotted contours) and when 
running is included (full contours). The area enclosed by the 
95% contour increases by around 20% when running is in- 
cluded. 



the variation in the constraints with scale is much less, 
as for instance is seen in Fig. [S] showing the ns-r plane. 
Indeed in this case we find that minimization of the area 
is not only unnecessary, but can actually be misleading, 
because parameters such as r can appear to be well con- 
strained even on scales where there is no meaningful data. 
The reason for this is that the restrictive class of mod- 
els under consideration force the spectra to behave in a 
particular way as they are extrapolated away from the re- 
gion where the bulk of the data lie, i.e. such constraints 
contain significant prior information as well as data infor- 
mation. This is also true to some extent for constraints 
on r in the running case studied earlier. 

Nevertheless, it is now interesting to compare the run- 
ning and no- running constraints. In the WMAP3 analy- 
sis the impression, from cornparison of the top- left panels 
of Figs. 12 and 14 of Ref. is of a huge deterioration 
in the constraints in the ns-r plane once running is in- 
cluded. The same is seen in Fig. 1 of Ref. [9] . However 
we now see that this is an artifact of the choice of scale 
where the constraints are portrayed. At the optimal scale 
there is some deterioration, due to parameter degeneracy, 
but the area increase within the 95% contour is only by 
about 20% as seen in Fig. [3 not by a factor of five as 
at fc = 0.002 Mpc~^. Consequently, inclusion of running 
leads only to a moderate deterioration in constraints on 
e and r\. 



V. INCLUDING MORE DATA 

We explore the robustness of our results by carrying 
out the same analysis for a broader compilation of data, 
now including shorter-scale CMB experiments and galaxy 
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FIG. 8: As Fig. 21 but now with the full dataset compilation. 
Note the modified axis ranges. 



correlation data from ACBAR [Tg, CBI [l|L VSA [ig. 
Boomerang [H, SDSS HH, and 2dFGRS T^j. 

Everything goes through as before. We find that the 
decorrelation scale of ns and running is 0.016 Mpc~^, 
which is not significantly different from WMAP3 alone. 
Though in general one would expect the decorrelation 
scale to change with dataset, in this case the WMAP3 
data are powerful enough that a shift is not seen. 

The constraints, particularly on r and hence e, do 
tighten significantly with the extra data, as is clear also 
in previous analyses including Ref. [ij . As an illustration 
of the results we obtained in this case, we show the ar- 
ray of constraints on the lowest-order e and r\ at different 
scales, Fig. [51 and the overlay of contours in the na,-r 
plane at the optimal scale, with and without running, in 
Fig.m 



VI. CONCLUSIONS 

We have investigated the issue of choice of scale in 
presenting marginalized parameter constraints. While 
we have focussed on WMAP constraints applied to in- 
flationary models, the same considerations apply much 
more widely. For example, in constraining density per- 
turbations using galaxy clusters, commonly the parame- 
ter (78, being the normalization of density perturbations 
smoothed on the scale 8/i~^ Mpc~"^, is quoted. However 
typically the normalization is best determined at a some- 
what larger scale than 8/i~^ Mpc, and marginalizing over 
parameters such as fio to quote constraints on erg can 
unnecessarily increase the statistical uncertainty on the 
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FIG. 9: As Fig. [71 but now with the full dataset compilation. 
Note the modified axis ranges. The area enclosed by the 95% 
contour increases by around 30% when running is included. 



for WMAP3 is 0.017 Mpc . This criterion can be used 
to define such a scale for any dataset compilation, and we 
found that the scale shifts hardly at all when other avail- 
able data are added to WMAP3. The optimal scale may 
also have some modest dependence on the choice of model 
parameters varied in a fit, for instance if non- negligible 
neutrino masses were included. One might even wonder 
whether it might be best to choose different scales for 
different observables, as the scalars and tensors are best 
constrained on quite different length scales, but we have 
not attempted this here. 

We have shown that the marginalized constraints on ng 
and r, or on e and ry, depend significantly on the choice of 
scale in the presence of running. By choosing the optimal 
scale, we find that constraints on those parameters are 
only mildly degraded by the inclusion of running as a pa- 
rameter, in contrast to the impression given if constraints 
are quoted at a non-optimal scale such as 0.002 Mpc""'^. 
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